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ABSTRACT 

We study the absorption along lines of sight toward high-z radio sources 
caused by the 21 cm transition of neutral hydrogen in the intergalactic medium 
(IGM) before reionization. Using semi-analytic methods, we compute the number 
density of observable features caused by both "minihalos" (bound objects that 
are unable to cool efficiently because of their small virial temperatures) and pro- 
togalactic disks. We show that both sets of features should be observable by the 
next generation of low-frequency radio telescopes, including the Low Frequency 
Array and the Square Kilometer Array, provided that sufficiently bright back- 
ground sources exist. The statistics of minihalo absorption features seen along 
lines of sight to radio-loud quasars offer a way to measure the evolution of the 
radiation background and the IGM temperature with cosmic time. Intersections 
with disks are much less common but cause significantly deeper absorption fea- 
tures that would be visible in the spectra of both radio-loud quasars and gamma- 
ray bursts (GRBs). The absorption feature caused by HI in the host galaxy of 
a GRB should be observable, offering a route to determine spectroscopically the 
burst redshift. 

Subject headings: cosmology: intergalactic medium, structure formation - galax- 
ies: active, radio, radio lines 



1. Introduction 

Despite its apparent simplicity, effective observational probes of the intergalactic medium 
(IGM) at high redshifts are difficult to find. Between recombination at z ~ 1000 and 
reionization at z ~ 6, the IGM was almost entirely neutral. In such a medium, the 
optical depth to Lya absorption at the frequency corresponding to a redshift z is r a « 
6.45 x 10 5 xhi[(1 + z)/10] 3//2 (Gunn & Peterson 1965), where xm is the neutral fraction and 
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where we have assumed the currently favored cosmological parameters (see below). This 
enormous optical depth renders the Lya forest, which is the most important tool for study- 
ing the (highly ionized) IGM at low and moderate redshifts, almost useless for detailed 
studies of the pre-reionization IGM. For example, the detection of a complete absorption 
trough shortward of the Lya resonance wavelength in the spectrum of the quasar SDSS 
1030+0524 may be the first evidence for the end of the reionization era (Becker et al. 2001). 
However, studies of rest-frame UV and optical quasar spectra are unable to place stronger 
constraints on the neutral fraction than x m > 10~ 2 (Fan et al. 2002). Attempts to constrain 
other parameters, such as the thermal state of the IGM, the radiation background, or the 
density structure, are also compromised by the extremely large optical depth at the Lya 
resonance. 

Clearly, we cannot obtain an in-depth understanding of the IGM before and during reion- 
ization through studies of the Lya forest. An alternative is to observe a weaker transition 
at the other extreme of optical depth: the 21 cm hyperfine transition of neutral hydrogen. 
To date, there have been several theoretical studies of emission or absorption of the neu- 
tral IGM against the cosmic microwave background (CMB) in this transition (Scott & Rees 
1990; Kumar et al. 1995; Madau et al. 1997; Tozzi et al. 2000; Shaver et al. 1999; Iliev et 
al. 2002). However, the predicted signals are extremely weak. Even the next generation of 
radio telescopes, such as the Square Kilometer Array 1 (SKA), will only be able to detect 
exceptionally large objects (e.g., Scott & Rees 1990; Kumar et al. 1995). Another strategy is 
to study the angular structure of the emission in order to constrain statistically the structure 
formation process at high redshift (e.g., Tozzi et al. 2000; Iliev et al. 2002). However, source 
confusion with low-z faint radio sources (Di Matteo et al. 2002) and free-free emission by 
the intergalactic medium (Loeb 1996) are likely to compromise such attempts. 

Alternatively, Carilli et al. (2002) have pointed out that luminous high-redshift radio- 
loud quasars serve as ideal background sources against which absorption by intervening gas 
can be seen: with a bright background source it becomes much easier to identify the weak 
absorption features expected to be produced by the IGM. Along such lines of sight one 
can map the "21 cm forest" of redshifted hyperfine absorption lines and hence study the 
neutral IGM over the redshift range 6 < z < 10. Using a simulation, Carilli et al. (2002) 
found that the high-z analog of the Lya forest of filaments and other overdense regions 
produces recognizable absorption features in the spectrum of a radio-loud quasar. They 
also demonstrated that such observations are feasible with an instrument similar to current 
designs for the SKA, provided that sufficiently bright radio-loud quasars exist at these early 
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epochs. 

Large scale numerical simulations, such as those studied by Carilli et al. (2002), are 
unable to resolve collapsed gas clouds on the smallest mass scales, although these clouds 
represent the most abundant clumps in the high-redshift IGM. In this paper we complement 
their study by exploring semi-analytically the statistics of these compact absorbing systems. 
Our semi-analytic approach allows us to examine the sensitivity of the results to changes in 
the input parameters and to identify the physical quantities that are best measured by radio 
absorption spectra. As we demonstrate in the following sections, 21 cm absorption spectra 
probe the radiation background and thermal state of the IGM at high redshifts and therefore 
could be instrumental in testing models of structure formation and reionization. 

The smallest bound objects in the IGM (commonly referred to as "minihalos") have 
virial temperatures below the threshold for atomic hydrogen line cooling, and so (in the ab- 
sence of a substantial abundance of molecular hydrogen, H 2 ) they cannot cool and collapse 
to form protogalaxies. Although minihalos cannot form stars, they are important for deter- 
mining the mean clumping factor of the IGM and for screening ionizing radiation from other 
objects (Barkana & Loeb 2002). Their properties are sensitive to the uncertain presence of 
molecular hydrogen; if H 2 forms in sufficient quantities it can act as a cooling channel for 
these minihalos (Barkana & Loeb 2001, and references therein), strongly suppressing their 
number density (and hence increasing the global star formation rate). Detection of these 
halos would provide important insights into the physics of hierarchical structure formation 
prior to reionization. 

Halos with virial temperatures above the hydrogen line cooling threshold will collapse to 
form stars. Thus we expect some fraction of all lines of sight to penetrate protogalactic disks. 
Such absorption systems would be the pre-reionization analogs of damped Lja absorbers. 
While disk intersections are rare, they can provide important information about the state 
of the neutral hydrogen in high- ,2 galaxies, including the distribution of disk masses, star 
formation rates, and gas temperature. 

In §2 we derive the optical depth to redshifted 21 cm radiation through both the dif- 
fuse IGM and collapsed objects. We show our numerical results in §3 and conclude in 
§4. Throughout the discussion we assume a flat, A-dominated cosmology, with density pa- 
rameters Q m = 0.3, Q\ = 0.7, Qb = 0.05, in matter, cosmological constant, and baryons, 
respectively. In the numerical calculations, we assume h = 0.7, where the Hubble constant 
is H = 100/i km s" 1 Mpc" 1 . 
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Model 



2.1. The Diffuse IGM 



The optical depth of the neutral IGM to redshifted hyperfine absorption is (Field 1959a; 
Madau et al. 1997) 
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Here /ip is Planck's constant, k is Boltzmann's constant, uq = 1420.4 MHz = vj{\ + 2) 
is the rest-frame hyperfine transition frequency, A 10 = 2.85 x 10 -15 s _1 is the spontaneous 
emission coefficient for the transition, T$ is the spin temperature of the IGM, Tcmb — 
2.73(1 + 2;) K is the CMB temperature at redshift z, and nni is the neutral hydrogen density 
at this redshift. In the second equality, we have assumed sufficiently high redshifts such 
that H(z) w H QU 2 (l + zfl 2 (which is well-satisfied in the era we study, z r > 6). Because 
reionization is expected to be sudden (e.g., Gnedin 2000) we take a neutral fraction xm = 1 
throughout our discussion. 

The HI spin temperature (i.e., the excitation temperature of the hyperfine transition) 
is (Field 1958) 

rp _ Tcmb + y a T a + y c T K 

s ~ i + y a + yc ' U 

Here Tk is the HI kinetic temperature and T a is the Lya color temperature (see Madau et 
al. 1997 for a precise definition). In the conditions pertaining to both the IGM and collapsed 
objects, T a = Tk (Field 1959b). The third term describes collisional excitation of Ts, while 
the second describes the Wouthuysen-Field effect, in which Lya pumping couples the spin 
temperature to the temperature of the radiation field (Wouthuysen 1952; Field 1958). The 
coupling from collisions is 

_ C10 
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where C±o is the collisional de-excitation rate of the (higher-energy) triplet hyperfine level 
and = hpVo/k = 0.068 K. We use a fit to the results of Allison & Dalgarno (1969) 
for C10 at Tk < 1000 K, with an extrapolation C\q oc T^ ' 33 to higher temperatures. The 
Wouthuysen-Field coupling constant is 
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where P\q is the indirect de-excitation rate of the triplet level due to absorption of a Lya 
photon followed by decay to the singlet level. For a diffuse Lya background, Madau et al. 
(1997) showed that 

P w « 1.3 x 10~ 12 J_2i s" 1 , (5) 

where J_2i is the intensity of the background radiation field at the Lya frequency in units 
of 10~ 21 erg cm" 2 s _1 Hz -1 sr _1 . Within the diffuse IGM, the gas density is small enough 
that collisions do not strongly couple T5 and Tk- However, collisions can effectively couple 
the two quantities within collapsed objects. Lya pumping couples the two temperatures in 
a density-independent manner, becoming effective when J_ 2 i > 1. Below we show results 
for several values of J-21, because the UV radiation field before reionization is an unknown 
quantity, likely varying strongly with both position and time. For reference, the simulation 
used by Carilli et al. (2002) yields J- 2 i(z = 10) « 0.8 and J- 2 i(z = 8) « 4. 

The spin temperature and optical depth therefore depend on the kinetic temperature 
of the IGM. Once Thomson scattering with the CMB becomes inefficient at the thermal 
decoupling redshift z d ~ 140, the IGM cools adiabatically as the universe expands (Barkana 
& Loeb 2001) until the first objects collapse. As soon as this occurs, radiative (and mechan- 
ical) feedback from these objects quickly heats the IGM. Because the radiation background 
prior to reionization is still unconstrained by observations, estimates of the temperature rely 
on simulations. We choose a particular form for the temperature evolution that roughly 
matches the mass- weighted temperature found in the simulation of Carilli et al. (2002): 

T 1GM (z) = l -°- 252+5 - 5 K. (6) 

This simple form has two convenient properties. First, heating begins at z ~ 18, close to the 
time when the first stars are expected to form. Second, it yields Tiqm(z = 6) = 10 4 K, close 
to the temperature of a photo-ionized medium at approximately the time when reionization 
may have occurred. We neglect any spatial variation in the temperature field, as is likely to 
occur in filaments and other overdense regions. For contrast, we include below some results 
for zero heating and for a temperature that is smaller than that given in equation (6) by a 
factor of 10. 



2.2. Minihalos 

In cold dark matter models, structure forms hierarchically (i.e., low-mass halos collapse 
first and the characteristic mass scale increases with cosmic time). The smallest halos cannot, 
however, form stars because they cannot cool efficiently: without metals, atomic hydrogen 
line cooling is only effective for halos with a virial temperature T vir > 10 4 K. Here, T vir is 
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defined as in equation (26) of Barkana & Loeb (2001). We use the term minihalo to refer to 
collapsed objects with T vir < 10 4 K. This condition sets the upper mass limit for minihalos; 
the lower mass limit is determined by the Jeans mass or, more precisely, the time-averaged 
Jeans mass (Gnedin & Hui 1998) using the IGM temperature given in equation (6). 
If molecular hydrogen is present, cooling becomes more effective (Ricotti et al. 2002), but 
most models predict that even a small UV background will efficiently dissociate H 2 in these 
halos (Haiman et al. 1997; Ciardi et al. 2000; Haiman et al. 2000). 

We compute the optical depth r u of the 21 cm transition along a line of sight through 
a minihalo with impact parameter a. For the remainder of this section, v refers to the 
frequency of the background radiation at the position of the minihalo (i.e., not including 
the cosmological redshift between the observer and the minihalo). Each collapsed object is 
surrounded by an accretion shock. Inside the shock, we assume that the gas overdensity 
equals that of the dark matter, which we take to have the universal NFW profile of Navarro, 
Frenk, & White (1997). In reality, gas pressure will flatten the core of the gas profile relative 
to the dark matter, but we find that lines of sight passing through the affected regions are 
so rare that they make only a negligible difference to our final results. 

Outside of the virialization shock, each halo is surrounded by a region of infalling gas. 
In an Einstein-de Sitter universe, the structure of the infall region can be described by a self- 
similar solution (Bertschinger 1985). At the redshifts we consider (6 < z < 10), deviations 
from an Einstein-de Sitter universe are small, and the solution is approximately valid. The 
Bertschinger (1985) profile fixes the shock to be at a constant overdensity relative to the 
background medium: the shock radius r s h « 4r v ; r /3, with a slight dependence on cosmology 
and redshift. The virial radius r vir is defined in equation (24) of Barkana & Loeb (2001). 
Similarly, the turnaround radius (i.e., the radius at which the gas is at rest with respect to 
halo) is r ta « 3.85r vir . We scale the velocity field to match the Hubble flow and the gas density 
to match the mean IGM density at large radii. Unfortunately, the Bertschinger (1985) profile 
has a much steeper central cusp (p oc r -2 25 ) than the NFW profile. We therefore continue 
to use an NFW profile within the accretion shock in order to avoid seriously overcounting 
the number of high-optical depth lines of sight. As a result, the density jump condition at 
the shock is only satisfied to ~ 10% for the halos we study. However, because all of the 
quantities of interest are integrated along the line of sight, this error makes little difference 
to our results. 

The HI atoms are assumed to have a Maxwellian distribution at each point. For r < r s h 
we set T K = T vir , while for r > r sh , we set T K = T IG m- Note that we therefore ignore 
temperature variations within the shocked region. The Bertschinger (1985) solution gives 
a postshock temperature T sh ~ T vir /3 and thus slightly increases the optical depth. We 
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also ignore the heating of the gas as it is compressed in the infall region; in the adiabatic 
limit, this leads to at most a factor ~ 2.5 error in the temperature. However, because the 
overdensity in the infall region is not large, the effects on the total optical depths are small. 
The spin temperature is then determined locally via equation (2). 

Once the halo profile is specified, the optical depth at a frequency v along a line of sight 
with impact parameter a from the halo center is 

_ 3h P c 3 A w [°° n HI {r) f [v(u) - v LOS (a, R)} 2 

Tv 32nkv 2 T s (r)^b(r) eXP \ b 2 (r) 

Here, r 2 = a 2 + R 2 , v{v) = c{y — z/ )/Vo, ^los^, R) is the infall velocity projected along the 
line of sight, and b 2 (r) = 2kTx(r)/m p is the Doppler parameter of the gas. For v sufficiently 
far from uq, the exponential factor is non-zero only in the diffuse IGM where Tk, and Ts 
are constants. At these frequencies, the integral becomes tihi/[TsH(z)] due to the Hubble 
flow, and we recover equation (1). 

Sample optical depth profiles for a = 5 x 10 6 M minihalo at z^ = 10 are shown 
in Figure 1. The solid curves assume a heated IGM and J_2i = while the dotted lines 
show analogous results for J_ 2 i = 10. Each set of profiles assumes a = 0.3, 1.3, and 3r vir , 
from top to bottom. In all cases, the optical depth at large v — u is simply that of the 
IGM, which is strongly suppressed for a strong T S -T K coupling through a Lya background. 
Careful examination of the a = 3r vir curves reveals that r„ < tigm for moderately large 
v — vq. This occurs because the velocity gradients within the infall region exceed that of the 
Hubble flow, and so for a range of frequencies the absorbing column is actually smaller than 
that of the diffuse IGM. 

As is evident from Figure 1, the line width varies with impact parameter. For a < 
r v ; r , the dominant contribution to t v comes from shocked gas within the halo, and once 
absorption due to the diffuse IGM is removed, the line profile is close to a Gaussian with 
a width determined by T v - lT . For a ^> r v ; r , the line of sight does not pass through shocked 
gas and so the profile is approximately a Gaussian of width Tigm- For a > r vir , the line 
of sight passes close to the shock, and the profile is best fit by an effective temperature 
T eff = 0.5(Tigm + T vir /3). [The factor of 1/3 results from the infall velocities at the shock in 
the Bertschinger (1985) profile; see the paragraph above equation (7).] The observed FWHM 
(defined relative to the diffuse IGM absorption) is then 

where T eS is chosen as described above and Zh is the redshift of the minihalo. Note that 
Figure 1 shows intrinsic line profiles (i.e., without the cosmological redshift between the 
observer and the halo). 
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Another quantity of interest is a measure of the total absorption due to the halo, or the 
"equivalent width." We define the intrinsic equivalent width through the relation 



int 



poo poo 

/ (1 - e~ Tv )dv - (1 - e~ nGM )dv, (9) 

Ju J U() 



where tiqm is given by equation (1). This measures the amount of flux absorbed by the 
minihalo over and above that absorbed by the diffuse IGM. Because of the cosmological 
redshift, the observed equivalent width (Az/) b s = (Az/) int /(l + zy L ). 

The observed equivalent width as a function of impact parameter is shown for several 
z h = 8 halos in Figure 2. Solid curves show results for M h = 10 7 , 10 6 , and 10 5 M from top 
to bottom , with J_ 2 i = and a heated IGM. The dotted line shows results for M h = 10 6 M 
and no IGM heating, while the dashed line shows results for the same mass, full IGM heating, 
and a radiation background J„ 2 i = 10. The steep drop at a/r vir ps 1.3 occurs because of 
the density drop past the accretion shock. Note that the IGM temperature has little effect 
on the equivalent width. Smaller halos have larger central equivalent widths (because such 
halos are cooler, and t v oc T^T^^ 2 ) but smaller equivalent widths at the most probable 
impact parameters because their column density is smaller. A cold IGM slightly increases 
the equivalent width at large a because more of the absorbing gas is located within the infall 
region of the halo, where the density is above the mean IGM density. On the other hand, the 
radiation background has a pronounced effect on the equivalent widths, particularly at large 
a. By effectively coupling T s to T vir or T IGM in the outskirts of the halo where collisions are 
ineffective, such a radiation field drastically decreases the optical depth contributed by these 
regions. 



2.3. Galactic Disks 

Halos with T vir > 10 4 K can cool through atomic line radiation and are expected to 
collapse and form protogalactic objects. We assume here that these objects form disks. We 
adopt the exponential disk structure of Mo et al. (1998), in which 

E < a > = ^wr""' (io) 

Rd = -^=jdAr vir . (11) 

Here £ is the surface density of the disk, R d is the disk scale length, M d is the disk mass, 
jd is the specific angular momentum of the baryons in units of their initial value upon halo 
collapse, and A is the spin parameter. Mo et al. (1998) found that jd = 1 was required in order 
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to match observational data (implying no substantial loss of baryonic angular momentum 
to the dark matter), given the fact that N-body simulations show that A has a lognormal 
distribution with a mean A = 0.05. We will assume for simplicity that all disks have this 
mean spin parameter A. 

The primary contribution to HI absorption in galaxies comes from the cold neutral 
medium. We assume that a fraction of the total baryonic mass of the halo, (fl^/fl^Mh, 
is in this phase. (The remainder can be in stars, the hot interstellar medium, or a hot 
component of the halo.) We assume for simplicity that the cold gas is smoothly distributed 
throughout the disk (although in reality it is likely to be bound up in small, dense clouds 
dispersed throughout the disk). We take a spin temperature Ts = 10 3 K, as inferred from 
observations of damped Lya absorbers at lower redshifts (Carilli et al. 1996). We further 
assume that the spin and kinetic temperatures are coupled, as expected for the radiation 
background in a typical galaxy. The optical depth is then 



T v = TiGM + T exp 
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In this expression, we have assumed a random distribution of disk inclinations and have taken 
the average path length through an inclined disk (so that So — > 7r 2 /4So). We have ignored the 
modifications to tigm caused by the velocity gradient and overdensity of the infall region, 
which differs slightly from the Hubble flow (and therefore, at a given velocity offset, the 
column of neutral hydrogen differs from that of the unperturbed IGM). Typically, the large 
optical depths of lines of sight through galaxies make such subtle differences unimportant. 
However, note that disks are likely to be in complicated environments, both because of 
feedback from the collapsed object and because high-redshift protogalaxies are expected to 
be highly clustered. Feedback can heat and ionize the region around the disk, reducing the 
IGM absorption well below the nominal tigm- The equivalent width can be calculated from 
equation (9); with our assumption of a constant T s = T K within each system, the relation 
becomes 
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3. 



Results 



3.1. 



Number Densities 



Using our model we can compute the density of lines of a given peak optical depth or 
equivalent width in a spectrum. For a quantity C (either r or equivalent width), the number 
of systems intersected with C > Co per redshift interval is 



Here dr/dz is the comoving length element per unit redshift, duh/dM^ is the halo mass 
function (or number of halos at redshift z with mass between and M^+dM^ per comoving 
volume), and A(Mh, z, Co) is the cross-sectional area of a halo (in physical coordinates) with 
the specified mass and redshift subtended by lines of sight for which C > Co- For minihalos, 
M min = Mai and M max is determined by the cooling criterion. For disks, M min is determined 
by cooling and M max = oo. In both cases, we use dn h /dM h given by the Press-Schechter 
formalism (Press & Schechter 1974), with the modifications suggested by Jenkins et al. 
(2001) in order to better match the mass function found by numerical simulations. For disk 
calculations, we again assume a random distribution of disk inclinations; the effective area 
of each disk is therefore 4/7r 2 times the area of the disk viewed face-on. In this case the area 
is simply 



Figure 3 shows dN(> tq) / dz, the number of objects intersected per redshift interval with 
a central optical depth greater than r . The dashed, solid, and dotted curves show results for 
z = 10, 8, and 6, respectively. All minihalo calculations assume a heated IGM and J_ 2 i = 1- 
The density of minihalo lines decreases with cosmic time because the minimum halo mass 
increases as Tigm increases. Note that the evolution with redshift is particularly pronounced 
for large r . 

It is clear that minihalo absorption systems are much more common than disks at low 
optical depths, but the rare high-optical depth features are dominated by collapsed disks. 
The flat disk distribution results from the exponential dependence of the optical depth on 
impact parameter. The primary uncertainty in our disk calculation is T s . A single cold cloud 
in the Milky Way typically has T s ~ 100 K (Dickey & Lockman 1990), but observations of 
moderate-redshift absorbers suggest Ts ~ 10 3 K (Carilli et al. 1996). If T$ is smaller than we 
assume, the optical depth through the disk increases substantially [r oc T s 3 ^ 2 ; see equation 
(13)]. However, as equation (16) shows, the distribution varies only logarithmically with Tq, 




(15) 




(16) 
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so the corresponding change in dN(> r )/dz will not be large. For example, in the range of 
optical depths shown in Figure 3, setting T s = 10 2 K increases the number density of lines 
by a factor of only ~ 2. 

Figure 4 shows the variation of dN(> r )/dz with IGM properties at z — 8. The solid 
curves assume J_2i = 1, the dotted curves assume J_ 2 i = 0, and the dashed curves assume 
J_2i = 10. Within each set, the left line assumes our standard IGM temperature and the 
right line assumes that Tigm is an order of magnitude below the value given by equation (6). 
The steep rise as r approaches tigm is caused by the divergence of the relevant cross-section 
for each halo in this limit. We see that the radiation background can be inferred from the 
absorbing continuum while Tigm can be estimated from the statistics of deep lines. At small 
r , the distribution is primarily determined by the radiation field J_2i- On the other hand, 
dN(> r )/dz depends much more strongly on T IG m than on J_ 2 i at large r . Such large 
optical depths require a line of sight passing close to the halo center, where the gas density 
is high enough for collisional coupling between Ts and Tk to be effective, even without a 
substantial radiation field. The temperature dependence arises because a cool IGM has a 
smaller Mm and thus a larger number of low- mass halos with small T vir (or large To). 

Carilli et al. (2002) examined lines of sight through a cosmological simulation and es- 
timated that dN(r > 0.02, z = 10) /dz ~ 50 and dN(r > 0.02, z = 8)/dz ~ 4. For their 
simulation, the average radiation background evolves between J_2i ~ 1-4 over this redshift 
range. Figure 4 shows that their z = 8 result lies in the range we would predict for this 
radiation field, but their z — 10 point is ~ 4 times higher than our prediction in Figure 
3. [If we choose Tigm to precisely match the mass-averaged temperature at z — 10 in the 
simulation of Carilli et al. (2002), the discrepancy is reduced to a factor of ~ 3.] A difference 
in results should be expected given that the two calculations are complementary and do not 
consider the same population of absorbers. The baryonic mass per particle in the simulation 
was 10 5 ' 7 M (Gnedin 2001), while even the largest of our minihalos have baryonic masses 
~ 10 6 7 M . Our formalism therefore probes objects of much smaller mass than does the 
simulation; on the other hand, the simulation probes the filamentary structure of the IGM 
which is not adequately described by our formalism. Furthermore, the simulation includes 
non-uniform heating and radiation. Overall, it appears that filaments and minihalos provide 
comparable contributions to the line statistics, with the relative importance of filaments 
increasing as redshift increases. 

Another useful quantity is the number of intersections per redshift interval with observed 
equivalent widths greater than a fixed value, dN(> EW)/dz. In Figure 5 we show this for 
the same redshifts and IGM characteristics as in Figure 3. Again, disks only dominate the 
number counts for very large equivalent widths. In contrast to the optical depth statistics, 
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we find little evolution of dN(> EW)/dz with redshift. This is caused in part because our 
definition of the equivalent width removes the effect of tigm (which decreases rapidly with 
cosmic time). In addition, the small, high optical depth halos at high-z yield narrow lines and 
therefore do not have substantial equivalent widths. As a result, the simultaneous increases 
in Mm and minihalo density with cosmic time appear to roughly balance each other. We 
do not show the dependence of the equivalent width statistics on J_ 2 i and Tigm because it 
exhibits the same trends as dN(> r )/dz, although at a weaker level. 

3.2. Synthetic Spectra 

Our model can also be used to generate Monte Carlo realizations of absorption spectra 
for illustrative purposes. The probability P(z) that a line of sight intersects an object in a 
redshift interval dz centered on z is given by equation (15), with A(Mh, z) simply detemined 
by the maximum impact parameter allowed in the calculation. (We find that choosing 
« max = 5r vir is sufficient to include all observable absorption lines from minihalos.) To 
generate a spectrum, we choose redshift intervals Az small enough so that P(z) < 0.01. 
For each bin, we generate a random number to determine if an intersection takes place. If 
one occurs, we choose randomly using the probability distribution duh/dM^ and then 
randomly select the impact parameter with an equal probability per cross-sectional area of 
the halo. Note that this method does not include clustering of minihalos and in principle 
allows overlap between minihalos centered in nearby bins (although in practice such overlap 
is highly insignificant). 

Sample transmission curves generated in this way and smoothed to 1 kHz resolution are 
shown in Figure 6. Here the transmission T is defined as T = 1 — exp(— r v ). In each panel, 
we show the portion of a spectrum from z = 8.25 to z — 7.75. Each panel is generated from 
the same set of intersected halos, but assumes a different radiation field in the calculation 
of the optical depths (J_2i = 0, 1, and 10, from top to bottom). The figure explicitly shows 
that by coupling T s and T K , a radiation field can strongly decrease both the absorption of 
the diffuse IGM and of the minihalos. [For reference, Tigm (z = 8) = 3162 K in our model.] 
The strongest absorption is caused by a 2 x 10 6 M minihalo at z — 8.2, with a = 0.35r v i r . 
Close inspection reveals that most absorption lines are surrounded by a small region with 
transmission greater than that of the diffuse IGM. Such features result from the steep velocity 
gradients in the infall region of a halo, which cause the optical depth to fall below tigm by 
small factor (see Figure 1). 

Similarly, one can generate realizations of spectra including disk intersections. Because 
these intersections are rare (see Figure 3), we performed a Monte Carlo procedure on the 



-13- 



interval between z = 8.25 and z = 7.75 sixteen times in order to find a line of sight passing 
through a disk with r > 0.1. [Figure 3 shows that for disks, dN(r > 0.1) /dz 0.09 at 
z = 8. On average, we would therefore expect one such intersection per 22 realizations of this 
redshift interval.] In fact, the resulting system has r = 3.76. We then added this system to 
an independently generated minihalo field, and we show the resulting spectrum in Figure 8. 
In order to contrast the effects of disk and minihalo intersections, we switch to logarithmic 
axes and plot the flux decrement, 1 — T = exp(— t v ), instead of the transmission. The strong 
feature is caused by the disk, while the much weaker features to the right are caused by 
minihalos. Although such lines of sight are obviously rare, we see that disks can in general 
be clearly distinguished from the minihalo population because they have much larger peak 
optical depths and equivalent widths. Note that the line width here is a direct result of our 
assumption that Tk = 10 3 K. Most damped Lya absorbers at lower redshifts show more 
than one velocity component of somewhat narrower intrinsic width (Briggs et al. 2001; Lane 
& Briggs 2001). 



4. Discussion 

We have shown that absorption at the 21 cm transition of neutral hydrogen in both 
minihalos and protogalactic disks at high-redshifts can cause non-negligible absorption in 
the low-frequency spectra of high-redshift radio sources. Using a semi-analytic approach, we 
have found that the density of minihalo absorption features is primarily a function of the 
kinetic temperature of the IGM (which determines the minimum mass of the minihalos) and 
the radiation background (which couples the spin and kinetic temperatures of the hydrogen 
gas). We predict that the average radiation background can be measured by the continuum 
absorption and the number density of very weak lines, while the IGM temperature can be 
measured by the density of narrow, relatively strong absorption lines. Optical depths reach 
ro ~ 0.02 for objects common enough to appear along typical lines of sight. Studies of 
absorption radio spectra therefore offer a probe of the era before reionization and can help 
constrain the history of early star formation. We find, unsurprisingly, that disks are only 
rarely intersected but tend to have much higher optical depths. 

The principal shortcoming of our semi-analytic approach is its inability to accurately 
describe the filamentary structure of the IGM and the spatial clustering of halos. Using a 
simulation, Carilli et al. (2002) have also studied 21 cm absorption along lines of sight to 
high-redshift radio sources. While the mass resolution of their simulation was too coarse 
to include the minihalo population, they were able to capture the geometry and structure 
of filaments. They find that both the typical optical depths and the number densities of 
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filament absorption lines are comparable to or slightly larger than our estimates for those 
from minihalos. An observed absorption spectrum will then be a mixture of lines from these 
two populations of absorbers. The clustering of minihalo lines may yield information on bias 
or feedback from nearby cooled objects, but a proper study of such an effect will require 
numerical simulations or more detailed semi-analytic models (e.g., Scannapieco & Barkana 
2002). 

Can these signals be detected by future radio telescopes? First, note that the relevant 
frequency range is 100 — 200 MHz, corresponding to z ~ 13-6; thus, a low-frequency in- 
strument is necessary. The transmission curves in §3.2 ignore noise. In a real observation, 
the noise level will be determined by a combination of the telescope characteristics and the 
brightness of the background source. Because the specifications of the next generation of 
telescopes are not yet settled, we have chosen not to present results specifically tailored to 
any particular instrument. Instead, we can simply estimate the required source brightness, 
S'min, in order to observe a single absorption feature, given the telescope's effective area A e g, 
system temperature T sys , and channel width Az/ C h, for a specified signal-to-noise ratio S/N 
and integration time t. In the small To limit and under the assumption that Az/ C h <C Az/ b s , 
the result is 

Current designs for the SKA call for A e s/T SJS ~2x 10 3 m 2 K 1 , while the Low Frequency 
Array 2 (LOFAR) is expected to have A cS /T sys ~ 4.5 x 10 2 m 2 K _1 . A channel width of 1 kHz 
should be achievable with LOFAR while maintaining a total spectral coverage ~ 4 MHz and 
with SKA while maintaining a total spectral coverage of ~ 10 MHz. 

Do objects of this brightness exist at sufficiently high redshifts? Haiman & Loeb (2001) 
have shown that, in the context of hierarchical cold dark matter models, black holes with 
masses > 10 8 M Q can plausibly exist even beyond z ~ 10. While we currently have no ob- 
servational constraints on this model, Carilli et al. (2002) have also presented a plausibility 
argument for the existence of luminous radio-loud quasars at z > 7. By extrapolating the ra- 
dio galaxy luminosity function to high redshifts using the observed decline of optical quasars, 
they estimate > 2000 sources on the sky of sufficient brightness to be useful. Furthermore, 
they point out that a bright radio galaxy at z — 5.2 has already been detected (van Breugel 
et al. 1999), and there is no obvious environmental reason for a high- ,2 cutoff. 

Because of their high optical depths, disks could be detected along lines of sight to much 
fainter sources. In this case, equation (17) must be modified to apply to large optical depths 



2 See http://www.lofar.org/index.html. 
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through the substitution r — > 1 — exp(— r ). For example, detecting a disk with To = 1 given 
the telescope parameters in equation (17) would require an average source flux ~ 0.4 mJy. 
Extreme starburst galaxies and faint active galactic nuclei could satisfy this criterion and of 
course will be more common than the bright quasars useful for minihalo studies. However, 
the disk intersection probability is sufficiently small that only a small fraction of all lines of 
sight will pass through intervening disks, so searches using these sources are unlikely to be 
practical given the long integration times required. Another particularly interesting type of 
source is a gamma-ray burst (GRB). In this case, the rarity of disks is not a limit because 
the host galaxy is necessarily present at the redshift of the burst. (While the host galaxy of 
a faint AGN or a starburst may also be visible, the gas is much more likely to be ionized 
on large scales around such sources.) The smooth power-law spectrum of GRBs is an ideal 
continuum for these observations. The typical peak flux of a GRB is < 1 mJy (Ciardi & 
Loeb 2000); however, the low-frequency spectra (< 10 GHz) of GRBs typically suffer from 
synchrotron self-absorption (Piran 2000), reducing the flux at 1420 MHz by a factor of a 
few from the peak value (e.g., Galama et al. 2000). GRBs are therefore too dim for the 
study of minihalo absorption features, but absorption by the host disk should be observable 
provided that r > 2. Such a detection could confirm the redshift of the burst and provide 
information on the neutral hydrogen content of its host galaxy. 
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Fig. 1. — Optical depth profiles for a = 5 x 10 6 M Q minihalo at Zh = 10 embedded 
in a heated IGM. Each set of profiles assumes a = 0.3, 1.3, and 3r vir , from top to bottom. 
The solid curves assume a heated IGM and J_ 2 i = 0, while the dotted lines show results for 

J -21 = 10. 
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Fig. 2. — Observed equivalent width as a function of impact parameter a. Solid curves show 
results for M h = 10 7 , 10 6 , and 10 5 M from top to bottom, with J_ 2 i = and a heated IGM. 
The dotted curve shows results for = 10 6 M and no IGM heating, while the dashed 
curve shows results for the same mass, IGM heating, and a radiation background J_ 2 i = 10. 
All curves assume Zh — 8. 
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Fig. 3. — Number of systems with observed central optical depth greater than To intersected 
per redshift interval. The dashed, solid, and dotted curves show results for z = 10, 8, and 
6, respectively. The two sets of lines describe minihalos and protogalactic disks. Minihalo 
calculations assume J_ 2 i = 1 and a heated IGM. 
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Fig. 4. — Number of systems with observed central optical depth greater than r intersected 
per redshift interval at z — 8. The solid curves assume J_ 2 i = 1, the dotted curves assume 
J_2i = 0, and the dashed curves assume J_2i = 10. Within each set, the leftmost curve 
assumes our standard IGM temperature and the rightmost curve assumes Tigm is an order 
of magnitude below that of equation (6). 
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Fig. 5. — Number of systems with observed equivalent width greater than a threshold in- 
tersected per redshift interval, as a function of the threshold value. The dashed, solid, and 
dotted curves show results for z = 10, 8, and 6, respectively. The two sets of lines de- 
scribe minihalos and protogalactic disks. All minihalo calculations assume a heated IGM 
and J_2i = 1. 
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Fig. 6. — Simulated transmission along a line of sight to a distant radio source assuming 
absorption by intervening minihalos. All panels are generated using the same realization of 
the minihalo mass field and assume a heated IGM. The panels show J_ 2 i = 0, 1, and 10, 
from top to bottom. 
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Fig. 7. — Simulated flux decrement profile along a sightline to a distant radio source chosen 
to intersect a protogalactic disk. The accompanying minihalo mass field assumes a heated 
IGM and J_ 2 i = 1. 



